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ABSTRACT 

We present chemical and radiative transfer models for the many far-IR 
ortho- and para-H^O lines that were observed from the Orion-KL region in 
high resolution Fabry-Perot (FP) mode by the Long Wavelength Spectrometer 
(LWS) on board the Infrared Space Observatory (ISO). The chemistry of the 
region was first studied by simulating the conditions in the different known 
components of Orion-KL: chemical models for a hot core, a plateau and a 
ridge were coupled with an accelerated A-iteration (ALI) radiative transfer 
model to predict H2O line fluxes and profiles. Our models include the first 45 
energy levels of ortho- and para-H^O. We find that lines arising from energy 
levels below 560 K were best reproduced by a gas of density 3xl0 5 cm~ 3 at 
a temperature of 70-90 K, expanding at a velocity of 30 km s _1 and with a 
H2O/H2 abundance ratio of the order of 2 - 3 x 10~ 5 , similar to the abundance 
derived by Cernicharo et al. (2006). However, the model that best reproduces 
the fluxes and profiles of H2O lines arising from energy levels above 560 K 
has a significantly higher H2O/H2 abundance, 1 - 5 x 10~ 4 , originating from 
gas of similar density, in the Plateau region, that has been heated to 300 
K, relaxing to 90-f00 K. We conclude that the observed water lines do not 
originate from high temperature shocks. 
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1 INTRODUCTION 



High velocity gas was first detected at the centre of the Orion-KL region as broad wings on 
'thermal' molecular lines in the millimeter range and as high velocity maser features in the 
22 GHz line of H 2 (Genzel et al. 1981). These high velocity motions may be caused by 
mass outflows from newly formed stars. Many theoretical studies of the Orion region (e.g. 
Draine & Roberge 1982; Chernoff et all. 1982; Neufeld & Melnick 1987) have concluded that 
the rich emission spectrum of thermally excited water vapor should play a significant role 
in cooling the gas. 

The widespread distribution of water vapour around IRc2 has been probed with maps 
at 183 GHz (Cernicharo et al. 1990, 1994, Cernicharo & Crovisier 2005), the first time that 
its abundance was estimated in the different large-scale components of Orion IRc2. Harwit 
et al. (1998) analysed 8 water lines observed with the ISO LWS FP, concluding that these 
lines arise from a molecular cloud subjected to a magnetohydrodynamic C-type shock. From 
their modelling, they derived an H 2 to H 2 abundance ratio of 5x 10 -4 . However, the in- 
terpretation of the lines observed in the ~ 80" LWS beam and the determination of the 
water abundance in the different Orion components remains a long standing problem, due 
to two main issues: the complexity of the dynamical and chemical processes that are taking 
place within the region encompassed by the LWS beam, with outflows and several different 
gas components, and the need for H 2 collisional rates appropriate for the temperatures 
prevailing in shocks. 

A total of 70 water lines were detected in the ISO-LWS far-IR spectral survey of Orion- 
KL (Lerate et al. 2006), with typically 70 km s" 1 FWHM. The line profiles range from 
predominantly P-Cygni at shorter wavelengths to predominantly pure emission at longer 
wavelengths. At the shorter wavelengths, the heliocentric velocities of absorption compo- 
nents appear to be centred at ~ - 15 km s" 1 , consistent with the results found in the 
ISO Short Wavelength Spectrometer (SWS) wavelength range, shortwards of 45 /im (van 
Dishoeck et al. 1998; Wright et al. 2000). However, the radial velocities of the pure emission 
lines of H 2 peak at around +30 km s _1 , whereas the velocity of the Orion-KL quiescent 
gas is +9 km s _1 . 
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In a previous analysis of the H 2 lines observed by the ISO LWS-FP, Cernicharo et al. 
(2006) modeled lines from the first 30 rotational levels of ortho- an para-H 2 0, concluding 
that the lines mainly arise from a gas flow expanding at 25-30 km s _1 , and inferred a gas 
temperature of approximately 80-100 K and a H 2 0/H 2 abundance ratio of 2-3 x 10 -5 . This 
derived abundance was an average over the LWS beam and they suggested that water could 
be locally more abundant if the emitting region included warmer components which interact 
with the ambient gas. 

In the current work, a technique to distinguish between the different components has 
been applied to the final calibrated high-resolution spectra of the ortho- and para-H 2 lines 
from the Orion-KL region. The methodology is similar to that used to model the ISO-LWS 
high resolution spectra of the Orion-KL CO lines (Lerate et al. 2008). Chemical models of 
the physically distinct components are coupled with a non-local radiative transfer model. 
The description of the models is structured to follow the different components of the region: 
the hot core, plateau and ridge. Many references can be found for a description of the KL 
region components - a complete overview is given in Stahler and Palla (2004). Both models 
(chemical and radiative transfer) are highly configurable and have been applied to a variety 
of emitting regions, including molecular gas in outflows (Benedettini et al. 2006). 



2 OBSERVATIONS AND DATA REDUCTION 

The ISO-LWS FP observations were carried out between September 1997 and April 1998. 
The dataset consists of 26 individual observations making up a total of 27.9 hours of ISO 
LWS observing time in L03 Fabry-Perot full spectral scan mode. The dataset also includes 
16 observations accumulated over 13.1 hours in L04 Fabry-Perot line scan mode and 1 
observation in the lower resolution L01 grating scan mode. The instrumental field of view 
for all L03 observations was centred either on a position offset by 10.5" from the BN object 
(which is at 05 h 35 m 14.12 s - 05° 22' 22.9" J2000), or on a position offset by 5.4" from IRc2 
(which is at 05 h 35 m 14.45 s - 05° 22' 30.0" J2000), while most of the L04 observations were 
centred on IRc2. The LWS beam had a diameter « 80" (Gry et al. 2003). Processing of the 
LWS FP data was carried out using the Offline Processing (OLP) pipeline and the LWS 
Interactive Analysis (LIA) package version 10. The basic calibration is fully described in the 
LWS handbook (Gry et al. 2003). Further processing, including dark current optimisation, 
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de-glitching and removal of the LWS grating profile was then carried out interactively using 
the LWS Interactive Analysis package version 10 (LIA10; Lim et al. 2002) and the ISO 
Spectral Analysis Package (ISAP; Sturm 1998). A detailed description of the observations 
and of the data reduction process is given by Lerate et al. (2006). 

3 THE CHEMICAL AND RADIATIVE TRANSFER MODELS 
3.1 The Chemical Models 

The chemical model used to simulate the Orion KL plateau and ridge is described by Viti 
et al. (2004) and is the same as used by Lerate et al. (2008) to model the Orion-KL lines 
observed by the ISO-LWS Fabry-Perot. The chemical network is taken from the UMIST 
database (Le Teuff et al. 2000). The model uses a two phase calculation in which gravita- 
tional collapse occurs in phase I, with gas-phase chemistry and sticking onto dust particles 
with subsequent processing (hydrogenation and conversion of a fraction of CO into methanol) 
also occurring. In phase II, where evaporation from grains also took place, we simulated the 
presence of a non-dissociative shock by an increase of temperature (from 200 to 2000K, 
depending on the model) at an age of ~ 1000 yr, which is the dynamical timescale of the 
main outflow observed in the KL region (Cernicharo et al. 2006). The duration of the high 
temperature shock is about 100 years, representing the temperature structure of a C-shock. 
The gas is then allowed to cool. This temperature profile was adopted from the calculation 
of Bergin et al. (1998) who studied the chemistry of H2O and O2 in postshock gas. The 
treatment of the temperature increase is considered to be linear with time. The modeling 
allows a wide range of parameters to be varied in order to investigate a range of condi- 
tions. As with the CO modelling of Lerate et al. (2008), the main parameters that were 
varied in this analysis were: (i) the initial and final densities (ii) the depletion efficiency, 
hereafter called the freez-out parameter (iii) the cloud size (iv) the maximum gas tempera- 
ture and (v) the interstellar radiation field (ISRF). Note that the H2 column density is not 
a free input parameter but is calculated self-consistently as a function of size and gas density. 

We based the parameter grid on descriptions of the Orion-KL components found in 
detailed analyses of the region (Blake et al. 1987, Genzel and Stutzki 1989, Cernicharo et 
al. 1994). The main parameters adopted from these references were the sizes and densities. 
However, we used our own analysis of the continuum emission and molecular rotational 
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Model 


Density (cm 3 ) 


T shock (Kj 


T 9 as (K) 


mco% 


size (//) 


N(H 2 )(cm- 2 ) 


PL1 


3xl0 5 


200 


80 


40 


30 


4xl0 22 


PL2 


3xl0 5 


300 


90 


60 


30 


4xl0 22 


PL3 


3xl0 5 


500 


90 


60 


30 


4xl0 22 


PL4 


3xl0 5 


500 


90 


80 


30 


4xl0 22 


PL5 


3xl0 5 


1000 


100 


60 


30 


4xl0 22 


PL6 


3xl0 5 


1000 


100 


80 


30 


4x 10 22 


PL7 


3xl0 5 


2000 


100 


60 


30 


4xl0 22 


PL8 


3xl0 5 


2000 


100 


80 


30 


4xl0 22 


PL9 


lxlO 6 


300 


90 


60 


30 


2xl0 23 


PL10 


lxlO 6 


300 


90 


80 


30 


2xl0 23 


PL11 


lxlO 6 


500 


90 


60 


30 


2xl0 23 


PL12 


lxlO 6 


500 


90 


80 


30 


2xl0 23 


PL13 


lxlO 6 


1000 


100 


60 


30 


2xl0 23 


PL14 


lxlO 6 


1000 


100 


80 


30 


2xl0 23 


PL15 


lxlO 6 


2000 


100 


60 


30 


2xl0 23 


PL16 


lxlO 6 


2000 


100 


80 


30 


2xl0 23 


RG1 


lxlO 4 


no shock 


70 


40 


15 


9xl0 20 


RG2 


lxlO 4 


no shock 


70 


20 


15 


9xl0 20 


RG3 


5xl0 4 


no shock 


80 


10 


15 


5xl0 21 


RG4 


5xl0 4 


no shock 


80 


20 


15 


5xl0 21 



Table 1. Plateau (PL) and Ridge (RG) models and their parameters: density, maximum temperature reached by the gas in 
the shock simulation, gas temperature after cooling, the percentage of mantle CO (mco) given by the freeze-out parameter at 
the end of Stage I of the chemical model, and the approximate H2 column density (assuming ~ 1.6 xlO 21 cm -2 per magnitude 
of visual extinction). 

diagrams (from Lerate et al. 2006) to indicate gas temperatures. Table [1] list the main 
parameter sets used to compute the grid of chemical models that were used for the present 
H2O modelling and for the CO modelling of Lerate et al. (2008). For completeness we also 
list the derived H 2 column densities for each model. 

The chemical model used to simulate the extended gas in Orion KL is based on the 
same model used for the shocked gas, with the exception that no high temperature shock 
is included as an input. The model simulates the gas chemistry evolution of gas expanding 
with velocities of ~ 25-30 km s _1 , which is heated up to ~ 100 K on a timescale of 1000 yr, 
relaxing to 80 K (see Figured]). A grid of conditions was also investigated for six extended- 
gas models (Table [2]). Their initial parameters were based on the published literature on 
the Orion-KL extended gas (Tielens and Hollenbach (1985), Draine and Roberge (1982), 
Menten et al. (1990)). 



3.2 The ortho-H 2 and para-H 2 radiative transfer models 

As described by Lerate et al. (2008), the chemical model produces abundances that are 
used as inputs to the radiative transfer model SMMOL (Rawlings & Yates, 2001), along 
with values for the dimensions, density, dust temperature and radiation field intensity. The 
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Model Density (cm" 3 ) T gaB (K) mco% expansion gas velocity (km s 1 ) 



El 3xl0 5 80 80 25 

E2 3xl0 5 90 70 30 

E3 5xl0 4 100 60 35 

E4 3xl0 5 80 50 25 

E5 3xl0 5 90 40 30 

E6 5xl0 4 100 30 35 



Table 2. List of main parameters modified to investigate the extended-gas chemical model of Orion KL. 

SMMOL code has accelerated A-iteration (ALI) that solves the radiative transfer problem 
in multi-level non-local conditions. It starts by calculating the level population assuming 
LTE and takes an adopted radiation field as the input continuum and then recalculates 
the total radiation field and level populations, repeating the process until convergence is 
achieved. The resulting emergent intensity distribution is then transformed to the LWS 
units of flux (W cm -2 /im -1 ) taking into account the different beam sizes (slightly different 
for each detector) and is convolved with an instrumental line profile corresponding to a 33 
km s _1 FWHM Lorentzian (Polehampton et al. 2007), in order to directly compare with 
the observations. The input parameters include molecular data such as the molecular mass, 
energy levels, radiative and collisional rates and also the dust size distribution and opacity. 
The parameters also include the gas and dust temperatures of the object being modelled. 
We estimated the dust temperature from the LWS grating observation of Orion-KL (Lerate 
et al. 2006), fitting a black body function of 70 K. In order to reproduce the observed 
continuum flux level, we adopted a radiation field equivalent to 10 4 Habings (1 Habing unit 
corresponds to the integrated flux in the wavelength range from 91.2 to 111.0 nm of 1.6 x 
10~ 3 erg s _1 cm -2 ; Habing 1968). The adopted ortho-para ratio was 3, which was found by 
Barber et al. (2006) to be appropriate when T>50 K. The microturbulent velocity was set 
to 5 km s" 1 and expansion velocities from 15-40 km s _1 were considered for the shocked gas 
in the Plateau models. We estimated an error of less than 30% for the fit to the continuum, 
being the maximum percentage deviation below the observed continuum for wavelengths up 
to 120 yum and the maximum percentage deviation above the observed continuum for longer 
wavelengths. 

The o-H 2 molecular data were taken from public molecular databases (Muller et al., 
2001; Schoier et al., 2005). Table |3] lists the input parameters for the radiative transfer model 
that were modified with respect to the CO modelling. The number of energy levels was set 
to 45 and the number of radiative transitions to 158. The number of collisional transitions 
was 990 and 10 collisional temperatures were investigated, from 20 K to 2000 K. The H 2 0- 
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INPUT PARAMETERS 


VALUE 


MOLECULAR. WEIGHT 


18.0 


NUMBER OF RADIATIVE TRANSITIONS 


158 


NUMBER OF COLLISION PARTNERS 


1 


NUMBER OF COLLISION TEMPERATURES 


10 


NUMBER OF COLLISIONAL TRANSITIONS 


990 


ISUM: THE STATISTICAL EQUILIBRIUM EQUATION FOR LEVEL ISUM IS 
REPLACED BY THE EQUATION OF CONSERVATION OF NUMBERS 


45 


NK: NUMBER OF ENERGY LEVELS INCLUDING CONTINUUM LEVELS 


45 



Table 3. List of main 0-H2O and P-H2O parameters adopted as inputs in the radiative transfer models 

H 2 collisional excitation rates were based on the (scaled) H 2 0-He calculations of Green et 
al. (1993). Note that since the calcuations for this paper were performed new collisional 
rates for H 2 0-H 2 have become available (Faure et al. 2007; Dubernet et al. 2009). Table 1 
of Faure et al. shows that for temperatures below 300K the critical densities with helium 
can be factors of 2-5 greater than for rates with ortho and para hydrogen. However as our 
densities are at least 3 orders of magnitude below the reported critical densities, the radiative 
pumping terms are as important, if not more so, than the collisional pumping terms. We 
would therefore expect our results to be qualitatively valid within the observational errors. 
However one would expect the use of He as a collisional partner to be incorrect when trying 
to model Herschel HIFI data, because the much better angular and spectral resolution of 
the observations may reveal pockets of much higher density; in such cases the use of the new 
collisional rates for water will be necessary. 

4 RESULTS AND DISCUSSION 

The H 2 line profiles and fluxes were reproduced by coupling our chemical models (see Fig- 
ure [1] and Figure [2] for the final adopted chemical models) with SMMOL radiative transfer 
models. A large number of models were investigated, with simulations of different shock 
temperatures, see Tableland Table [2j In the models where the H 2 lines were reasonably 
well reproduced, an extended grid of conditions was then investigated to refine the line fits. 
An example of this is shown for the Extended gas model, whose investigated scenarios are 
summarised in Table IU The main parameters that were modified to refine the models were 
the Phase II gas temperature and the freeze-out parameter, which determines the amount 
of water frozen onto the grains at the end of phase I of the chemical model. 
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Figure 1. Time evolution H2O, OH and CO abundances corresponding to the extended gas chemical model E2 in Table [2] 




Figure 2. Time evolution of H2O, OH and CO abundances for the PL2 chemical model (see Table[T]for a description of the 
main model parameters). 
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Model 


Density (cm 3 ) 


T 9 as (K) 


gas velocity (km s 1 ) 


mH 2 (%) 


H2O/H2 abundance 


Wl 


3xl0 5 


80 


25 


70 


2.1 X 10~ 5 


W2 


3xl0 5 


90 


30 


60 


2.8 X 10~ 5 


W3 


5xl0 4 


100 


35 


50 


4.8 X 10~ 5 


W4 


3xl0 5 


80 


25 


10 


7.5 x 10~ 5 


W5 


3xl0 5 


90 


30 


80 


9.8 x 10" 5 


W6 


5xl0 4 


100 


35 


20 


4.8 x 10" 6 



Table 4. List of main parameters and resulting H2O/H2 abundances from the chemical models adopted to investigate the 
extended-gas component in Orion-KL. 
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Table 5. H2O line fluxes observed and predicted by the modelling. [1] Observed LSR velocity. Where two values are listed, these correspond to separate emission peaks. 
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Line profile fits from our radiative transfer modelling are shown in Figure |3] and Figure 
HI while Table |5] compares the observed and predicted line fluxes. The final column of Table 
[5] lists the x 2 values, defined as x 2 = J2iLi(fl ux - m °d(i) — f lux -obs(i)) 2 /N (bins), for the 
comparison between the modeled (flux_mod) and observed (flux_obs) line profiles. The line 
profile fits are excellent for the pure emission lines, while the lines with P Cygni profiles are 
reasonably well reproduced. 

Overall, our results can be summarised as follows: 

• Lines arising from energy levels below ps 560 K are best reproduced by 70-90 K gas of 
density 3xl0 5 cm -3 , expanding at approximately 25-30 km s" 1 ([3]). A graphical representa- 
tion of the time evolution of H 2 0, OH and CO is shown in Figured] corresponding to Model 
E2 in Table |2j Model W2 from Table [4] gives very similar results. 

• Lines from higher energy levels, both ortho and para, are best reproduced by warmer 
gas of density 3xl0 5 cm -3 , expanding at approximately 25-30 km s" 1 , which is heated up to 
300 K and then relaxes to 90-100 K (Figure H]). This corresponds to Model PL2 in Table [U 
which also reproduced the observed CO transitions with J up ^18 (Lerate et al. 2008). Note 
that this is simply the best x 2 fit: models where the shock temperature is higher, such as 
PL3, also provide a good fit to the line profiles. However models where the shock temperature 
is of the order of 1000 K, the derived water abundance did not reproduce the line profile. 
We should also emphasize that the parameter that most affected the fits to the water lines 
was the freeze-out rate, as this directly affects the fractional abundance of water (since the 
higher the freeze-out rate the more oxygen is hydrogenated as water on the surfaces of the 
grains before evaporating). Based on the investigated models, shock temperatures in the 
300-500 K range were found to produce acceptable fits to the line fluxes and profiles. 

• The PL2 model which fits the higher-excitation lines shown in Figure H] makes a negli- 
gible contribution to the lower-excitation lines shown in Figure [3J while the E2 model that 
fits the lower excitation lines makes a negligible contribution to the higher-excitation lines 
shown in Figure HI 

• The H 2 to H 2 abundance ratio is of the order of 2-3 x 10 -5 for model W2 and 1-5 x 
10~ 4 for model PL2. 

• For the 14 para-H 2 lines that show pure emission, the ratio of the predicted to observed 
emission flux is found to be 0.89 ± 0.39, while for the 10 ortho-H 2 lines that are purely in 
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Figure 3. Predicted line profiles from the PL2 radiative transfer model (solid lines), over-plotted on the observed ortho and 
para-H20 lines (dotted lines). The PL2 model corresponds to gas of density 3xl0 5 cm -3 , heated up to 300 K and relaxing to 
90-100 K, expanding at a velocity of 30 km s~ 1 (Tablc[TJ. The ordinate corresponds to the continuum-subtracted flux and the 
abscissa is the radial velocity in km s _1 . 



emission this ratio is also 0.89 ±0.39. We interpret this as observational evidence in support 
of the 3:1 ortho:para ratio that was adopted for the modeling. 

The results from our modeling of the observed far-IR ortho and para-H 2 lines from 
Orion-KL show that two different chemical models are needed to reproduce the H2O lines in 
the LWS wavelength range. For lines arising from energy levels below ~ 560 K, our modelling 
results are in agreement with the findings of Cernicharo et al. (2006); our chemical model of 
an expanding gas with velocity of 30 km s _1 at 70-90 K is able to reproduce the line fluxes 
and line profiles of both the ortho- and para-H^O lines, with an H2O/H2 abundance ratio of 
the order of 2-3 x 10~ 5 . For the water lines that exhibit P Cygni profiles, our current profile 
fits appear to provide an improvement when compared to the fits presented by Cernicharo 
et al. (2006). 



However, for lines arising from higher energy levels (above 560 K) the model that best 
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Figure 4. Predicted line profiles from the E2 radiative transfer model (solid lines) over-plotted on the observed ortho and 
para-H20 lines (dotted lines). The E2 model corresponds to gas of density 3xl0 5 cm -3 at 70-90 K expanding at a velocity 
of 30 km s _1 (Table [2j. The ordinate corresponds to the continuum-subtracted flux and the abscissa is the radial velocity in 
km s~ 1 . 



reproduces both the H2O line fluxes and their profiles is of warmer gas which is initially 
heated up to ~300 K and then relaxes to 90-100 K. The corresponding H 2 0/H 2 abundance 
ratio depends on the time stage within the PL2 model but is of the order of 1-5 x 10~ 4 , 
within reach of the value of 5 x 10~ 4 derived by Harwit et al. (1998) from their modelling 
of eight Orion-KL LWS-FP water lines. For the seven water lines in common, Harwit et 
al.'s mean observed/predicted line flux ratio was 1.9±1.5, versus a mean value of 1.4±0.4 
from our modelling. Our PL2 model was also found to reproduce well the observed LWS-FP 
CO transitions having J up ^18, interpreted as arising from the Plateau region within the 
extended warm component emission (Lerate et al. 2008). Note, however, that, for higher-J 
CO lines our work shows that a higher temperature gas is needed, in agreement with other 
authors (e.g. Watson et al. 1985), confirming the findings of Lerate et al. (2008) that the 
observed molecular emission arises from multiple components that differ in density and tern- 




perature. 

To conclude, we find that, taken together, Plateau region model PL2 and extended-gas 
region model E2 can match all of the measured far-IR water lines from Orion-KL. The 
main difference between the two modelled zones is that the Plateau region has warmer 
temperatures, with a consequent impact on the chemical evolution. As noted by Cernicharo 
et al. (2006), radiative pumping due to the strong IR dust continuum radiation field is 
sufficient to populate the higher excitation rotational water lines. The H 2 0/H 2 ratio in the 
extended-gas region is found to be ~ (2 — 3) x 10~ 5 , similar to the value found by Cernicharo 
et al. (1996) from their line modelling, but we find a significantly higher ratio, ~ (1—5) x 10~ 4 , 
in our Plateau region models that fit the profiles and fluxes of the higher-excitation water 
lines. 

Our present results should be taken together with those of Lerate et al. (2008), who 
analysed ISO LWS-FP observations of multiple rotational lines of CO and found that in 
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order to explain the emission from all of the the CO transitions, hot cores as well as shocked 
regions (with temperatures ranging from 300 to 1000 K) had to be present within the ISO 
beam. 
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